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Constraining the symmetry energy content of nuclear matter 
from nuclear masses: a covariance analysis 
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Elements of nuclear symmetry energy evaluated from different energy density functionals 
parametrized by fitting selective bulk properties of few representative nuclei are seen to vary widely. 
Those obtained from experimental data on nuclear masses across the periodic table, however, show 
that they are better constrained. A possible direction in reconciling this paradox may be gleaned 
from comparison of results obtained from use of the binding energies in the fitting protocol within 
a microscopic model with two sets of nuclei, one a representative standard set and another where 
very highly asymmetric nuclei are additionally included. A covariance analysis reveals that the 
additional fitting protocol reduces the uncertainties in the nuclear symmetry energy coefficient, its 
slope parameter as well as the neutron-skin thickness in 208 Pb nucleus by ~ 50%. The central values 
of these entities are also seen to be slightly reduced. 

PACS numbers: 21.30.Fe,21.65.Ef,21.60. Jz 


I. INTRODUCTION 

Constraining quantitatively the symmetry energy C® 
of nuclear matter at saturation density po and its den¬ 
sity slope L 0 (= 3po dC Qp P ^ Ipo) h as been a major focus of 
attention in present-day nuclear physics. They are fun¬ 
damental parameters in influencing the binding energies 
and stability of atomic nuclei and in shaping the isospin 
distributions there giving rise to neutron-skins (A r np ) in 
nuclei of large neutron excess [J-@|, where A r np is the 
difference between the rms radii for density distributions 
of the neutrons and protons in a nucleus. They are also of 
seminal importance in astronomical context. The inter¬ 
play of gravitation with the pressure of neutron matter 
P„(= 3poLo) at po is a determining factor in the radii of 
neutron stars Q. The dynamical evolution of the core¬ 
collapse of a massive star and the associated explosive 
nucleosynthesis depend sensitively on the slope param¬ 
eter To B @l- ^ a I S0 controls the nature and stability 
of phases within a neutron star, its critical composition, 
thickness and frequencies of crustal vibration Q and de¬ 
termines the feasibility of direct Urea cooling processes 
within its interior BIS. 

The value of the symmetry energy coefficient C 3 has 
been known for some decades to be in the range of 
~ 32.0 ± 4.0 MeV p~oMT3| . In a controlled finite-range 
droplet model (FRDM) from a fit of the observed nuclear 
masses, the uncertainty in its value shrunk appreciably, 
it was found to be = 32.5 ± 0.5 MeV S- Studying 
meticulously the double differences of symmetry energies 
estimated from experimental nuclear masses, Jiang et al 
|T~il find its value to be 32.10±0.31 MeV. Their procedure 
takes advantage of the fact that other effects like pairing 
and shell corrections in symmetry energy are well can¬ 
celled out from double differences of neighboring nuclei. 

Initial explorations on the symmetry slope Lo, how¬ 
ever, show wide variations. The suggested lower limit is 
~ 20 MeV, on the high side, it is ~ 120 MeV 0 ], Astro- 


physical observations of neutron star radii and masses 
provide 0 ] a value of Lq between 43 and 52 MeV, 
though. On the other hand, correlation systematics of 
A r np with nuclear isospin 0], isoscaling [0], nucleon 
emission ratios [3 , isospin diffusion 0 in heavy-ion col¬ 
lision, analysis of giant dipole resonance of 208 Pb [20l [2li| . 
giant quadrupole resonance in 208 Pb jH}, pigmy dipole 
resonance in 68 Ni and 132 Sn [2f| or of nuclear ground 
state properties using the standard Skyrme Hartree-Fock 
approach J24[, all yield values of Lq that differ substan¬ 
tially from one another [25j . Efforts have also been made 
to constrain the value of Lo by fitting the binding ener¬ 
gies of large number of deformed nuclei within the density 
dependent point coupling model (Hj, 0 and by consid¬ 
ering the pseudo data on the equation of state for the 
asymmetric nuclear matter or pure neutron matter in 
the fitting protocol of the Skyrme models 0, 0 as well 
as the density dependent meson exchange models j3Cj. 
These studies indicate that L 0 is in the range of 40 - 70 
MeV. 

Analysis of data on the nuclear masses in macroscopic 
models seems to contain the fluctuations in the value 
of Lq somewhat better. In a macroscopic nuclear model 
such as liquid drop, surface tension tends to favor a nearly 
uniform drop of neutrons and protons. For a nucleus, say, 
with large neutron excess it is energetically advantageous 
to distribute some neutrons out at the surface of reduced 
density (where the symmetry energy is lower) giving rise 
to a neutron skin. Imprints of symmetry energy coeffi¬ 
cient C 3 and its density slope Lq are thus encoded in the 
precisely known nuclear masses. Fitting the few thou¬ 
sand observed nuclear masses within the FRDM [3], the 
value of Lq was found to be 70±15 MeV. The surface 
symmetry coefficient C s , along with C 3 might also be 
used to estimate the value of L 0 10. Consistent with 
the values of and C s (as determined from double dif¬ 
ferences of symmetry energies estimated from experimen¬ 
tal nuclear masses 0] to be 58.91±1.08 MeV), a recent 
effort aided by microscopic calculations on the A r np of 
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TABLE I: The best fit values for the parameters of model-I and model-II. m a is the mass of a meson given in units of MeV. 
The masses of uj and p mesons are kept fixed to m„= 782.5 MeV and m p = 763 MeV and nucleon mass is taken to be M= 
939 MeV. Uncorrelated errors on the fitted parameters (upper line) and total correlated errors on them (lower line) are given 
inside the parenthesis for both the models. 


Name 

9° 


9p 

«3 


V2p 

Co 

m a 

model-I 

-10.62457 

13.8585 

12.077 

1.46285 

-0.9673 

28.33 

5.2056 

496.0067 


(0.00013) 

(0.0002) 

(0.045) 

(0.00024) 

(0.0024) 

(0.47) 

(0.0027) 

(0.0062) 


(0.246) 

(0.662) 

(2.60) 

(0.275) 

(3.66) 

(29.9) 

(3.21) 

(12.2) 

model-II 

-10.62123 

13.85989 

12.436 

1.46223 

-0.8566 

32.50 

5.3220 

495.8146 


(0.00011) 

(0.00017) 

(0.046) 

(0.00021) 

(0.0022) 

(0.49) 

(0.0027) 

(0.0054) 


(0.149) 

(0.262) 

(1.54) 

(0.290) 

(1.53) 

(18.1) 

(0.099) 

(8.23) 


heavy nuclei gives L 0 = 59 ± 13 MeV. This is 

not far from the one ( Lq = 65.5 ± 13.5 MeV) obtained 
from the nuclear equation of state (EoS) constructed [dij 
from some well-accepted empirical macroscopic nuclear 
constraints. Exploiting nuclear masses thus seems more 
or less to constrain better the symmetry elements C 3 and 
Lq of nuclear matter. 



FIG. 1: (Color online) The binding energy difference 
A B{X,Y) = B(X) - B(Y) for four different pairs of iso¬ 
topes are plotted against neutron-skin thickness A r np in the 
208 Pb nucleus for 16 different RMF models (See text for de¬ 
tails). The values of Pearson’s correlation coefficients r are 
also displayed. 

Energy density functionals (EDF) in microscopic mean 
field models, parametrized to reproduce the binding en¬ 
ergies of nuclei alongwith some other specific nuclear ob¬ 
servables do not, however, display such constraints on 
the symmetry elements of nuclear matter. The EoS con¬ 
structed out of the different EDFs show wide variations in 
the value of C 3 and Lq j35j, HQ]. Questions then arise how 
the information content of symmetry energy gets blurred 


in the exploration of nuclear masses in microscopic mod¬ 
els. For example in the study [37] of the binding energy 
difference A B of 132 Sn and 100 Sn (which is sensitive to 
the symmetry energy content) vis a vis the neutron skin 
A r np of 132 Sn for different sets of Skyrme EDFs, no no¬ 
ticeable correlation between A B and A r np was found. 
The different EDFs reproduce the binding energy differ¬ 
ence fairly well, but the value of A r np differ in a wide 
range. 

This quizzical finding leaves open the understanding 
on the incompatibility of the tighter constraints on the 
symmetry energy coefficients of nuclear matter extracted 
from experimental nuclear masses with the larger uncer¬ 
tainties on the values of C 3 and Lq as obtained from 
parametrizations of the microscopic mean field models. 
The purpose of this communication is to give a general 
direction in shading light on this seeming paradox. In 
an exploratory calculation, we show that the apparent 
irreconciliation can be eased if some highly asymmetric 
even-even spherical nuclei are additionally included in the 
fitting protocol of the optimization of the EDFs in mi¬ 
croscopic mean field models. The relativistic mean-field 
model (RMF) is chosen as the vehicle for the realization 
of our goal. 

To keep the matter simple yet pointed, we exploit the 
binding energy difference between four pairs of nuclei in 
different RMF models showing increasing asymmetry ef¬ 
fects: ( 68 Ni - 56 Ni), ( 132 Sn - 100 Sn), ( 24 0 - 16 O) and 
( 30 Ne — 18 Ne). The neutron rich 68 Ni and 132 Sn nuclei 
have asymmetries <5=0.176 and 0.242 respectively (<5 is 
the isospin asymmetry parameter (N — Z)/A); 24 0 and 
30 Ne have <5=0.3. The Ni and Sn isotopes are doubly 
closed shell nuclei. So also the O-nuclei, 24 0 is recently 
seen to be an unexpectedly stable doubly magic nucleus 
[H, . The Ne-nuclei have their neutron shells closed 

but have valence protons. The binding energy differ¬ 
ence between the two Ne-nuclei is expected to cancel the 
pairing and the possible core-polarization effects arising 
from the two valence protons partially. In Fig[l]the bind¬ 
ing energy difference between the four pairs of nuclei are 
plotted against A r np of 208 Pb, the A r np and the bind¬ 
ing energies being calculated for seven models of BSR 
family [40, [4l[ , NL3 [Hf , FSU Q and for seven models 
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of Density Dependent Meson Exchange (DDME) family 
[43j . The correlation coefficient for the Ni-pair is seen to 
be only 0.012, for the Sn-pair, it has increased to 0.586. 
For the O and Ne pairs, they are quite high, 0.980 and 
0.978, respectively. One can not but fail to notice the 
increasingly high correlation with increasing asymmetry, 
particularly for the latter two cases. 

The occurrence of strong correlation for the case of 
O and Ne pairs probably suggest that selective combi¬ 
nation of suitable binding energies of nuclei of low and 
high isospin may be ideally suited to better constrain 
the isovector part of the nuclear interaction. A covari¬ 
ance analysis uncovers correlations among the physical 
observables; revelation of this interdependence may be 
of great importance in assigning values to the physical 
observables of interest with statistical uncertainties. We 
aim to do that presently, the observables we explore are 
those that carry the imprint of isovector content of nu¬ 
clear interaction i.e. the symmetry energy (7° and the 
symmetry energy slope Lq along with the A r np of 208 Pb. 


II. EFFECTIVE LAGRANGIAN 


The effective Lagrangian density for the RMF model 
employed in the present work is similar to that of the FSU 
one [S|, I44t - l46| . The interaction part of the Lagrangian 
can be written as, 


£int — ^P 


1 €■ 

+ ^gpT-Pf, + + t z) A p 




/■v4 


^3 2 3 

~m 9 ° rn '’ a 

9lmluj^p v p v . 


1 


24M 2 ' 5 ' m ' (j4 + 24 C off'M") 2 


(1) 


where, Nd and N p are the number of experimental data 
points and the number of fitted parameters, respectively. 
A Oi is the adopted error and O exp and Of h (p) are the 
experimental and the corresponding theoretical values for 
a given observable. The values of A Oi are chosen in such 
a way that Xo(p) I s unity. The experimental observables, 
in our case, include only the values of the binding energies 
and charge radii of a selected set of nuclei. There may 
be values of parameters that are also likely, they can 
give a reasonable fit to the data, the reasonableness of 
the domain of parameter space being defined as y 2 (p) < 
Xo(p) + 1 |47H49| . A reasonable domain of parameters 
yields a reasonable set of results for each set of predicted 
observables. 

To explore the results of Fig|T]more critically, two mod¬ 
els (model-I and model-II) corresponding to different sets 
of fit-data are constructed. In model-I the binding ener¬ 
gies and charge radii of some standard set of nuclei ( 16 0 , 
40 Ca, 48 Ca, 56 Ni, 68 Ni, 90 Zr, 100 Sn, 116 Sn, 132 Sn, 144 Sm 
and 208 Pb) spanning the entire periodic table are taken 
as fit-data. In model-II, we have the same set of experi¬ 
mental observables, but with the addition of the binding 
energy difference A B of ( 24 0 , 16 O) and of ( 30 Ne , 18 Ne). 
The parameters of model-I and model-II are obtained by 
optimizing the objective function y 2 (p) as given by 
Eq. ©. 

Once the optimized parameter set is obtained the cor¬ 
relation coefficient between two quantities A and B , 
which may be a parameter as well as an observable, can 
be evaluated within the covariance analysis as, 


AAAB 

CAB = — - - - ; - 

V AA 2 AB 2 


(3) 


where, covariance between A and B is expressed as, 


It contains the usual Yukawa coupling between the nu¬ 
cleonic field (ip) and different meson fields, namely, 
isoscalar-scalar a (coupling constant g a ), isoscalar-vector 
u> (g u j), isovector-vector p {g p ) and photon Affi The pa¬ 
rameters K 3 , K 4 and (0 determine the strength of self¬ 
couplings of a and uj mesons. The self-couplings enable 
one to produce softening in the EoS of symmetric nuclear 
matter if needed. The ui-p cross-coupling with coupling- 
constant r] 2 p controls the density dependence of symme¬ 
try energy. 


III. FIT DATA AND MODEL PARAMETERS 


The values of the parameters entering the EDF of the 
RMF model are obtained from an optimal y 2 fit of the ex¬ 
perimental observables with the theoretically calculated 
values. The optimal values of the parameters of the EDF 
give the minimum y 2 (p) °f the objective function y 2 (p) 
(p referring to the parameter space), defined as, 


x 2 (p) = 


1 


Vb / 0 e»P _ Q th ( p ) 

N d -N p y{ A Oi 


( 2 ) 



Here, C a p is an element of inverted curvature matrix 
given by, 

<W = I(fAEl) . (5) 

2 V dp a dpp >po 

where, po represents the optimized set of parameters. 
The square of the error, AA 2 in A can be computed us¬ 
ing Eq. fi} by substituting B = A. Curvature matrix 
can also facilitate to locate the reasonable domain of pa¬ 
rameters 00 [p. The parameters for model-I and II 
corresponding to minimum value of the objective func¬ 
tion x 2 (p) (=Xo(p)) along with their uncorrelated and 
correlated statistical errors are listed in Table Q] The cor¬ 
related and uncorrelated errors are computed with and 
without contributions from the non-diagonal elements of 
the matrix C as given by Eq. [5j respectively. The cor¬ 
related errors on the parameters for the case of model- 
II are smaller than those obtained for the model-I indi¬ 
cating that the inclusion of the fit data on the binding 
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TABLE II: Observables O of different nuclei, adopted errors on them AO, their experimental values and the ones obtained 
for model-I and II. B and r c h refers to binding energy and charge radius of a nucleus respectively, and A B is binding energy 
difference of two isotopes of a nucleus as indicated. B and A B are in units of MeV and r c h in fm. 


Nucleus 

O 

AO 

Expt. 

model-I 

model-II 

16 0 

B 

4.0 

127.62 

127.781i0.990 

127.783i0.576 


r c h 

0.04 

2.701 

2.700A0.017 

2.699i0.013 

16 o, 24 o 

A B 

2.0 

41.34 

- 

40.995il.046 

ls Ne, 30 Ne 

A B 

2.0 

79.147 

- 

79.149il.296 

40 Ca 

B 

3.0 

342.051 

342.929il.064 

342.927i0.927 


T ch 

0.02 

3.478 

3.457i0.013 

3.455i0.010 

£- 

00 

0 

p 

B 

1.0 

415.99 

414.883i0.720 

414.751i0.541 


Teh 

0.04 

3.479 

3.439i0.007 

3.439i0.006 

56 Ni 

B 

5.0 

483.99 

483.752i2.495 

483.619il.646 


Teh 

0.18 

3.750 

3.695i0.025 

3.693i0.020 

68 Ni 

B 

1.0 

590.43 

592.294i0.784 

592.162i0.736 

90 Zr 

B 

1.0 

783.893 

782.855i4.833 

782.776il.621 


r c h 

0.02 

4.269 

4.267i0.009 

4.267i0.034 

100 Sn 

B 

2.0 

825.8 

827.987il.753 

827.757il.534 

116 Sn 

B 

2.0 

988.32 

987.169i0.946 

987.072i0.760 


Teh 

0.18 

4.626 

4.623i0.009 

4.623i0.008 

132 Sn 

B 

1.0 

1102.9 

1102.851il.146 

1102.631i0.856 


Teh 

0.02 

4.71 

4.711i0.011 

4.712i0.010 

144 Sm 

B 

2.0 

1195.74 

1195.834il.240 

1195.736il.287 


Teh 

0.02 

4.96 

4.956i0.009 

4.956i0.009 

208 Pb 

B 

1.0 

1636.446 

1636.457i4.301 

1636.383i0.917 


Teh 

0.02 

5.504 

5.530i0.012 

5.531i0.010 


energy differences constrain the model parameters bet¬ 
ter. In particular, the errors on the parameters g p and 
r] 2 p , which govern the isovector part of the effective La- 
grangian, are smaller for the model-II. The large error 
on the parameters K 3 and K 4 for both the models may 
be due to the fact that the fit data does not include any 
observable which could constrain the value of the nuclear 
matter incompressibility coefficient [46]. In Table HD dif¬ 
ferent observables Oi , adopted errors on them A Oi, their 
experimental values along with the results obtained for 
model-I and model-II using the corresponding best fit pa¬ 
rameters are listed. The values of Oi and A Oi, except for 
the A B of ( 24 0 , 16 O) and ( 30 Ne , 18 Ne) and r c h of 132 Sn, 
are exactly same as used in Ref. [52| . The experimen¬ 
tal data for A B of ( 24 0 , 16 O) and ( 30 Ne , 18 Ne) are taken 
from [53j and that for the r c h of 132 Sn from [HU ]. 


IV. RESULTS AND DISCUSSIONS 

We now compare the results obtained for the model- 
I and model-II to see upto what extent the inclusion of 
the experimental data on the binding energy differences 
between the pair of O and Ne nuclei can constrain the 
iso-vector part of the effective Lagrangian. In Fig. Ha) 
the reasonable domain for the parameters g p and r/ 2 p is 


TABLE III: The values for the binding energy per nucleon 
E/A, incompressibility coefficient K, Dirac effective mass of 
nucleon m*/m, symmetry energy coefficient (7° and density 
slope parameter of symmetry energy Lo for the nuclear matter 
evaluated at saturation density po along with the correlated 
errors on them obtained within the covariance analysis for the 
models I and II. The results for neutron-skin thickness A r np 
in 48 Ca, 132 Sn and 208 Pb are also presented. 


Observable 

model-I 

model-II 

E/A (MeV) 

-16.036 i 0.070 

-16.036 i 0.051 

K (MeV) 

210.12 i 27.87 

209.64 i 28.52 

po (fm -3 ) 

0.150 i 0.003 

0.150 i 0.003 

m* /to 

0.585 i 0.012 

0.585 i 0.010 

C° v (MeV) 

32.03 i 3.08 

31.69 i 1.51 

Lo (MeV) 

57.62 i 17.08 

55.63 i 7.00 

Ar np ( 48 Ca) (fm) 

0.191 i 0.036 

0.187 i 0.016 

A r np ( 132 Sn) (fm) 

0.266 i 0.070 

0.257 i 0.031 

Ar np ( 208 Pb) (fm) 

0.201 i 0.065 

0.193 i 0.030 


displayed. For this reasonable domain of the parame¬ 
ters the values of the symmetry energy slope parameter 
Lq and the A r np in the 208 Pb nucleus are displayed in 
Fig. Hb). The inclined and elongated shapes of the ellip- 
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FIG. 2: (Color online) The covariance ellipsoids for the 
parameters g p - r/ 2 P (upper panel) and the corresponding 
To - A r„ p (lower panel) for the model I (color blue) and model 
II (color red). The area inside the ellipsoids indicate the rea¬ 
sonable domain of the parameters. 

soids indicate that the correlations amongst g p - r] 2 P and 
Lq - A r np are strong. In fact, the values of the correla¬ 
tion coefficients (Eq. J3])) for these pairs of quantities for 
both the models turn out to be ~0.95. It is evident that 
the ellipsoids depicting the results for the model-II (color 
red) are narrower in comparison to those for the model-I 
(color blue). This is suggestive of the fact that the in¬ 
clusion of the binding energies for the 24 O and 30 Ne put 
tighter constraints on the isovector part of the effective 
Lagrangian density. 

Nuclear matter properties for model-I and model-II are 
compared in Table urn Errors on the entities describing 
the isoscalar behavior of nuclear matter ( E/A , K, pq and 
m*/m ) are pretty much the same for both the models 
concerned. For model-II, however, a significant improve¬ 
ment (by a factor ~ 2) on the spread of parameters like 
and Lq, which describe the symmetry behavior of 


nuclear matter, is achieved over model-I. Strikingly, the 
errors on C 3 and Lq for the model-II agree very well with 
the ones obtained for the SAMi Skyrme force j2f| which 
includes the variational EoS for the pure neutron matter 
as pseudo data in the fitting protocol. We also provide 
the values of A r np for 48 Ca, 132 Sn and 208 Pb nuclei in 
Table Mil The reduction in the errors on A r np for the 
model-II in comparison to those for the model-I are in 
harmony with the results depicted in Fig. [2] 

To this end, it may be mentioned that the correlations 
of neutron-skin thickness in 208 Pb nucleus with Lq has 
been studied within the covariance approach for different 
nuclear energy density functionals [551 ] . These correla¬ 
tions exhibit some degree of model dependence. In view 
of this, the present investigation should be extended to 
other energy density functionals. The value of Lq might 
also be constrained by including in the fitting protocol 
the experimental data for the iso-vector giant dipole po¬ 
larizability for the 208 Pb nucleus (48|. However, it must 
be reminded that the Lq is well correlated with the prod¬ 
uct of dipole polarizability and the symmetry energy co¬ 
efficient at the saturation density C p rather than the 
dipole polarizability alone [2lj . 


V. SUMMARY AND CONCLUSION 

To summarize, we have tried to give a critical look on 
the question why the information content on the elements 
of symmetry energy do not get transported properly in 
the EDFs when they are parametrized to reproduce the 
binding energies of representative nuclei in a mean-field 
model. The asymmetry parameters are then less con¬ 
strained. Working in the confines of RMF model, we 
find that better constraint can be achieved in the isovec¬ 
tor part of the effective Lagrangian density and hence in 
the said symmetry elements if one includes, in the fit¬ 
ting protocol of the microscopic model binding energies 
of some highly asymmetric nuclei. In this preliminary 
investigation, comparison of results with two sets of nu¬ 
clei, one a standard set and another with two additional 
nuclei of very high asymmetry (6 = 0.3) shows a pointer 
in this direction. The covariance analysis yields the re¬ 
duction in the uncertainties on the values of symmetry 
energy and its slope parameter by ~ 50% in comparison 
to those obtained without the inclusion of the additional 
nuclei in the fit data. Similar effects are also observed on 
the values of A r np in 48 Ca, 132 Sn and 208 Pb nuclei. A 
more detailed analysis along this direction is underway. 
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